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Abstract 

The BeppoSAX Catalog has been very recently pub- 
lished. In this paper we analyze - using the Maximum 
Likelihood (ML) method - the duration distribution of 
the 1003 GRBs listed in the catalog with duration. The 
ML method can identify the long and the intermediate 
duration groups. The short population of the bursts is 
identified only at a 96% significance level. MC simu- 
lation has been also applied and gives a similar signif- 
icance level; 95%. However, the existence of the short 
bursts is not a scientific question after the Compton 
Gamma-Ray Observatory's observation. Our minor re- 
sult is this well-known fact that in the BeppoSAX data 
the short bursts are underrepresented, mainly caused 
by the different triggering system. Our major result 
is the identification of the intermediate group in the 
BeppoSAX data. Therefore, four different satellites 
(CGRO, Swift, RHESSI and BeppoSAX) observed the 
intermediate type Gamma- Ray Burst. 
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1 Introduction 



(jKouveliotou et al. 1993T) identified two types of Gam- 
ma-Ray Bursts (GRBs) based on durations, for which 
the value of Tg (the time during which 90% of the 
flucnce is a ccumulated) is smaller o r lar ger than 2 s, re- 
spectively. jMukheriee et al. 19981) and lHorvath 1998l )| 
made a suggestion about the third type of GRBs, 
which have intermediate (2-9 s) du ration. Later several 
papers confirmed this suggestion ( Hakkila et al. 2000l : 
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Balasteeui et al. 200ltlRaianiemi fc Mahonen 20021: IHorvath 2002|j 



Hakkila et al. 20031 : IHorvath et al. 20061 : IChattopadhvav et al. 2007t ).| 

All of these works were base d on BATSE observa tions. 
Lately, analysis of Swif t ( Horvath et al. 20081) and 



RHESSI (|Ripa et al. 2007tlRipa et al. 20091 ) data have 
found a similar group structu re (for RHESSI the \ 2 
method found only two groups (jRipa et al. 20071 ). how- 
ever, the Maximum Likelihood (ML) method re veals a 
significant intermediate group (jRipa et al. 20091 )). 

Since The Gamma-Ray Burst catalog obtained with 
the Ga mma Ray Burst Monito r (GRBM) aboard Bep- 
poSAX (jFrontera et al. 2009af ) has been published re- 
cently, this catalog is worth to be analyzed, whether 
the statistical methods can show us similar or different 
groups of GRBs. 

In Sect. 2 we analyze the duration distribution of the 
GRBs observed by BeppoSAX. In Sect. 3 we discuss the 
instrumcntational details of our analysis. In Sect. 4 we 
present the conclusions. 



2 Analysis of the duration distribution 



In the BeppoSAX catalog (jFrontera et al. 2009al : iFrontera et al. 2009b 
there are 1082 GRBs, of which 1003 have duration in- 
formation (see Figure 1. for the distribution). 

2.1 Maximum Likelihood calculations 



Similarly to (jHorvath 20021 ) we fit the logTgo distri- 
bution using Maximum Likelihood (ML) method with 
a superposition of k log-normal components, each of 
them having 2 unknown parameters to be fitted with 
N = 1003 measured points in our case. Our goal is to 
find the minimum value of k suitable to fit the observed 
distribution. Assuming a weighted superposition of k 
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log-normal distributions one has to maximize the fol- 
lowing likelihood function: 




Fig. 1 The duration distribution of the BeppoSAX 1003 
bursts. 
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where wi is a weight, /; a log-normal function with 
log Ti mean and 07 standard deviation having the form 
of 
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and due to a normalization condition 
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We used a simple CH — h code to find the maximum of 
Lk- Assuming only one log-normal component the fit 
gives L\ max = 5951.895 but in the case of fc=2 one gets 
L2max = 6011.355. See Table 1. for the best parame- 
ters. 



In (jHorvath et al. 2008f) we summarize the ML 
method and refer to ( Kendall fc Stuart 197fj| ) the con- 
fidence region of the estimated parameters, which is 
given by the following formula, where L max is the max- 
imum value of the likelihood function and Lq is the 
likelihood function at the true value of the parameters: 
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Based on this equation we can infer whether the ad- 
dition of a further log-normal component is necessary 
to significantly improve the fit. We make the null hy- 
pothesis that we have already reached the true value 
of fc. Adding a new component, i.e. moving from fc 
to fc + 1, the ML solution of L^ max has changed to 
L(k+i)max - but Lq remained the same. In the meantime 
we increased the number of parameters with 3 (wk+i, 



Table 1 The best parameters for the two log-normal fit of 
the GRB (observed by BeppoSAX) duration distribution. 



Fig. 2 The duration distribution of the BeppoSAX bursts 
with the Gaussian functions (long bursts with dashed line). 



Duration(logTgo) &(logTgo) 



shorter 0.62 0.62 306.2 

long 1.45 0.40 696.8 
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logTk+i and <J(k+i))- Applying Eq. (j4|) on both 



and L 



(k-\-l)max 



we get after subtraction 



2(L ( 



k-\-l)max 







(5) 



For fc = 1 Li2max is greater than L\ max by more than 
59, which gives for ;\J an extremely low probability. It 
means the two log-normal fit is a really better approxi- 
mation for the duration distribution of GRBs than one 
log-normal. The longer duration group (with a centroid 
about 28 seconds) has 69% of the burst population and 
the shorter one (with a centroid about 4.2 seconds) has 
31%. 

Thirdly, a three-log-normal fit was made combin- 
ing three fk functions with eight independent param- 
eters (three means, three standard deviations and two 
weights). The highest value of the logarithm of the like- 
lihood {L^ max ) is 6015.585. For two log-normal func- 
tions the maximum was L,2max = 6011.355. The maxi- 
mum thus improved by 4.23. Twice of this is 8.46 which 
gives us the probability of 3.7% for the difference be- 
tween L 2ma x and L 3max is being only by chance. 

2.2 Monte Carlo simulations 

One can check the 0.037 probability, using a Monte- 
Carlo (MC) simulation. Take the two-log-normal distri- 
bution with the best fitted parameters of the observed 
data, and generate 1003 numbers for Tqq whose distri- 
bution follows the best fit of the two-log-normal distri- 
bution. Find the best likelihood with five free parame- 
ters, two means, two dispersions and two weights (only 
one is independent). Make a fit with a three-log- normal 
distribution (eight free parameters, three means, three 
dispersions and two independent weights). Take the 
difference between the two logarithms of the maximum 
likelihoods and that gives one number in our MC sim- 
ulation. 

This procedure have been carried out for 1000 sim- 
ulations. There were 49 cases when the log-likelihood 
difference was more than the one obtained for the Bep- 
poSax data (4.23). Therefore, the MC simulations do 
not confirm exactly the numerical result obtaining by 
Eq. ((5]), however, it gives a similar (4.9%) probability 
for the statement that a third group is only a statisti- 
cal fluctuation ( for more details see (jD'Agostini 2003 ; 
D'Agostini 20041 )). 

The conclusion is the same; the third group popula- 
tion is not confirmed in high significance level. There- 
fore, there is a chance the third log-normal is not 
needed. This third population is the shortest dura- 
tion GRB group. However we already know from 
the BATSE sample the short duration group is exist. 



Therefore from the SAX data, using only the duration 
information we can confirm the existence of the short 
bursts with only 95-96% significance level. 



3 Discussion 

There is no question that the short bursts exist, since 
many satellites identified short (less than two seconds 
long) and spectrally harder bursts. In this paper we 
are using only the duration information of the bursts. 
The ML method was not able to identify all the three 
subgroups in the BcppoSax data at a high signifi- 
cance level (once again, using only Tqq). The third 
gro up in our analysis w as the shortest in duration. 
As (jFrontera et al. 2009af ) pointed out short bursts are 
less pronounced and displaced toward higher durations 
with respect to BATSE. This discrepancy is due to the 
lower efficiency of the GRBM trigger system to short 
GRBs, which, for almost the entire BeppoSax mission 
duration, used 1 s as short integ r ation time. More de 



tails can be seen in (Band 2003; Guidorzi et al. 2001 



Frontera et al. 19971 ). 

However, in the duration distribution two compo- 
nents were found at a very high significance level. Sur- 
prisingly, analyzing the duration distribution observed 
by BeppoSax GRBM one can find the long and the in- 
termediate duration population. The short population 
can be seen only with low significance level. 



4 Conclusions 

In the BATSE data many scientific groups identified 
three types of bursts; short hard , long soft and interme- 
diate duration very soft bursts (IM ukhcri ee et al. 1998 : 

Horvath 19981: lHakkila et al 200d iBalastegui et al. 200ll:l 



Raianiemi fc Mahonen 20021 iHorvath 2002tlHakkila et al 2003l :| 

Horvath et al. 200d IChattopadhvav et al. 20071). In 
the Swift faorvath et al. 20081) and RHESSI (|Ripa et al. 2007l :| 
Ri'pa et al. 20091 ) data recent works identified the same 
group structures. 

In this paper we analyzed a fourth data set - ob- 
served by the BeppoSAX satellite. The duration dis- 
tribution cannot be well fitted with one Gaussian com- 
ponent. The two Gaussian fit is much better and the 
three Gaussian fit is better only at a 95-96% signifi- 
cance level. However, we surely know from previous 
works that short (third group in our analysis) bursts 
exist. 

The existence of the intermediate type burst is still 
not widely accepted. In this paper we showed that in 
order to fit the BeppoSAX data duration distribution a 
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second component of intermediate duration is needed. 
However, the physical existence of the intermediate 
GRBs can be still questionable. ( Hakkila et al. 20031 ) 
argued that the statistical existence of the intermediate 
group in the BATSE sample is caused by instrumental 
effects. Instrumental biases can play a relevant role 
in shaping the observed duration distribution and the 
different fractions of short GRBs observed in the vari- 
ous catalogues reflects this. In the case of the GRBM, 
an important bias is given by the relatively long (1 s) 
short integration time adopted through out most of the 
mission lifetime (IFrontera et al. 2009al ). Due to this, 
the claim for a physically separate class of intermedi- 
ate GRBs requires a word of caution and the upcoming 
catalogues of new missions, such as Fermi/GBM, will 
help to clarify this issue. 
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